Driven Rydberg atoms reveal quartic level repulsion 
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The dynamics of Rydberg states of a hydrogen atom subject simultaneously to uniform static 
electric field and two microwave fields with commensurate frequencies is considered in the range of 
small fields amplitudes. In the certain range of the parameters of the system the classical secular 
motion of the electronic ellipse reveals chaotic behavior. Quantum mechanically, when the fine 
structure of the atom is taken into account, the energy level statistics obey predictions appropriate 
for the symplectic Gaussian random matrix ensemble. 
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Quantum chaos considers correlations between the 
quantal properties of dynamical system and its classi- 
cal behavior. In particular level statistics of quantum 
systems chaotic in the classical limit should generically 
obey Random Matrix Theory (RMT) predictions H ac- 
cording to the famous conjecture of Bohigas, Giannoni 
and Schmit The link between RMT and chaotic sys- 
tems has been quite fruitful for quantum chaos studies 
(see for reviews ||[§)- 

Depending on the symmetries of a given strongly 
chaotic system, statistical properties of its quantum spec- 
trum fall into one of the three classes known from RMT: 
orthogonal, unitary and symplectic. The orthogonal class 
is typically associated with Hamiltonians invariant with 
respect to some generalized time-reversal symmetry (re- 
ferred to also as an antiunitary symmetry), the corre- 
sponding statistical ensemble of random matrices is re- 
ferred to as Gaussian Orthogonal Ensemble (GOE). In 
the absence of any such a symmetry the corresponding 
class is known as Gaussian Unitary Ensemble (GUE). Fi- 
nally energy levels of half- integer spin systems with no ge- 
ometrical symmetries and obeying the generalized time- 
reversal symmetry T squaring to minus unity (T 2 = — 1) 
reveal twofold Kramers degeneracy and their statisti- 
cal properties pertain to Gaussian Symplectic Ensemble 
(GSE). 

The conjecture has been tested on a number of theo- 
retical models (we refer the reader to reviews J3|,|| rather 
than numerous original papers). Experiments in atomic 
systems have been restricted to the orthogonal univer- 
sality class Q only. To break antiunitary symmetries in 
atomic species one needs, e.g., a static magnetic field, 
nonhomogeneous on the atomic scale || which is quite 
hard to realize. For that reason GUE type statistics have 
been observed experimentally for microwave billiards (so- 
called wave chaos experiments) where the time-reversal 
symmetry could be broken by applying some ferrite de- 
vises ifl. 



As shown by us recently atomic systems also al- 
low to generate GUE type statistics provided one uses 
not only the static fields but also microwaves. Then, in- 
stead of considering properties of eigenvalues of a given 
Hamiltonian H (which is not possible for time-dependent 
microwave perturbation) one considers quasienergies of 
the Floquet operator H = H — ihd/dt well defined for a 
periodic driving JlO[ . For appropriately chosen combina- 
tion of elliptically polarized microwaves and a static elec- 
tric field, the (generalized) time-reversal symmetries are 
broken and GUE type statistics may be observed. Even 
for weak fields when Floquet states may be thought of 
as the perturbed principal quantum number n hydro- 
genic manifold, the typical splitting between levels may 
be of the order of a few MHz making the experiment fea- 
sible. A RMT type level statistics is a manifestation of 
the classically chaotic behavior as observed (via classi- 
cal perturbation theory) for the motion of the electronic 
ellipse (the so-called secular motion). 

With GOE or GUE type of statistics realized for 
atomic species it seems natural to ask whether it is pos- 
sible to observe also statistics corresponding to the third 
universality class - the symplectic ensemble. The aim of 
this letter is to provide an example of such a situation. 
Such a behavior has not yet been observed, as far as we 
know, neither for atomic nor for billiard type systems. 

To observe the effect one needs an interaction explic- 
itly depending on the electron half-integer spin, an inter- 
action sufficiently intense to relatively strongly perturb 
the levels. For hydrogen atoms two obvious candidates 
are the Zeeman effect and the spin-orbit coupling. In 
the former case, however, there is no time-reversal sym- 
metry squaring to minus unity. The latter seems quite 
weak in the semiclassical Rydberg states regime. For that 
reason traditionally the spin-orbit coupling is neglected 
while considering strongly externally perturbed Rydberg 
states. 

It is another story if we restrict ourselves to small ex- 
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ternal perturbations and consider the combined effect of 
these perturbations and the spin-orbit coupling. Our pre- 
vious experience |||9| still indicates that the secular mo- 
tion of the electronic ellipse may be strongly chaotic even 
for weak (but judiciously chosen) combination of the ex- 
ternal fields. Then, if also spin-orbit coupling is taken 
into account, one may hope to observe GSE type statis- 
tics. 

Thus we consider the very same model as before, i.e., 
the hydrogen Rydberg atom driven by microwaves and 
placed in a static uniform electric field and we add the 
spin-orbit interaction. The configuration of the external 
fields must be chosen so that: the time-reversal symmetry 
with T 2 = — 1 is preserved, any other geometrical sym- 
metry is broken and the dynamics is irregular classically. 
To fulfill the first condition we have to use microwaves 
of linear polarization, in contrast to our earlier studies 
HU ■ The Hamiltonian of the system in atomic units is 
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where a is the fine structure constant. H\ is the external 
perturbation 



Hi = E ■ r + Fx cos ut + F'z cos 2wi 
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gives the lowest order relativistic corrections |lE . Among 
them the most important for the present contribution is 
the spin-orbit coupling. In the above formulae E denotes 
a static electric field vector, F and F' are amplitudes 
of two microwave fields with frequencies u> and 2u> and 
polarized along the Ox and Oz axes, respectively. 

The Hamiltonian (H) is time reversal- invariant, with 
time-reversal operator squaring to minus unity. Thus 
eigenvalues of the Floquet Hamiltonian must reveal a 
twofold Kramers degeneracy [pj. The role of F' field is 
to break the geometrical symmetry which is preserved by 
(|l|) if F' = 0. Indeed if we choose, for F' = 0, the Oz axis 
such that a static field component parallel to it vanishes, 
E z = 0, then the parity operator Tl z = exp (iirS z )P z 
(where P z is a reflection with respect to Oxy plane) with 
the property H 2 — — 1 commutes with H. This uni- 
tary symmetry splits the quasienergy spectrum into two 
blocks with identical eigenvalues (Kramers degeneracy). 
Each block is in itself a hermitian matrix, i.e. its eigen- 
values are expected to obey GUE but not GSE statistics 
if the underlying classical motion is chaotic || . 

Since both the microwaves and the static electric field 
are assumed to be weak, the model may be analyzed per- 
turbatively both classically and quantum mechanically. 
The latter is straightforward in the effective Hamilto- 



nonvanishing order (first order for the static electric field 
due to a linear Stark effect as well as for the spin-orbit 
interaction and second order in microwave amplitudes F 
and F'). The resulting matrices of dimension 2n§ with 
no being the principal quantum number of the perturbed 
manifold studied have been diagonalized with standard 
routines. We refer the reader to for practical details 
described there for a similar model. 
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FIG. 1. Poincare surface of section (at $ = 0) of 
the classical secular motion, Eq. (^), for the hydrogen 
atom in static electric and microwave fields for the energy 
(Hq S + 1/2) /Eq = — 1 and with the amplitude and frequency 
of the microwave field Fq/Eo = 2, uio = 1.304, respectively. 
The coordinates used for the plot are the scaled angular mo- 
mentum Lo — L/no and its canonically conjugate angle ^. 
Panel (a) corresponds to the orientation of the static electric 
field vector 9 = 7r/2, tp — 7r/4 while panel (b) is related to 
the system additionally perturbed by F' field with the am- 
plitude Fq 2 /Eo = 5 (the static field orientation is 8 = 0.317T, 
tp — 7r/4). Note that, for the parameters chosen, not the 
whole (Lo,^) space is accessible. 

In an analogous way one may construct an effective 
classical Hamiltonian averaging over the phase of the field 
and the fast motion along the electronic ellipse and con- 
sidering the slow secular motion of the ellipse itself only. 
For the classical analysis we neglect the relativistic cor- 
rection part Hi (the spin-orbit interaction and Darwin 
term, <5(r), do not have direct classical analogs). For 
the purpose of the present analysis we assume the mi- 
crowaves to be off resonant thus the perturbation calcu- 
lations follow the Lie approach of |p^|,p^ ,p| rather than 
that appropriate for resonant driving [BL 

The resulting classical effective Hamiltonian (first or- 
der in E and second order in F and F') takes a form 

#o° ff = -~ + #i c ?e>o, F 2 , F ' 2 , S ; L , *, M , $) (4) 



where Hq S 



n 2 H cS 



, F 
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nian formalism |12| with the fields taken in the lowest 



ngf, E = UqE, 
luq = tiqUj, Lq — L/riQ and Mq = M/tiq are scaled vari- 
ables (L, M being the angular momentum and its pro- 
jection on the Oz axis while ^ and $ are angles con- 
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jugate to L and M, respectively). Observe that clas- 
sical dynamics depends only on the reduced energy via 
scaled variables (no dependence on no). Fig. [l]a shows 
an example of Poincare surface of section obtained for 



with respect to -l/2ng" between -2 • 10" 6 and 3 ■ 1CT 7 
while for the broken symmetry case between —3 ■ 10~ 6 
and -1CT 6 . 



(Hf + l/2)/E 



-1 with Fq = 0, 



Fi/Eo = 2, 



ujq = 1.304 and orientation of the static field vector 
9 = 7r/2, <p = tt/4 (where 9, ip are usual spherical angles). 
In Fig. |l|b there is a similar plot but for the case with bro- 
ken the parity symmetry, i.e. for the same parameters as 
previously but with Fq /Eq = 5 and 9 — 0.317T. Clearly 
the motion is predominantly chaotic in both cases. 
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FIG. 2. Nearest neighbor spacing distribution and spec- 
tral rigidity, A3 statistics, for hydrogen atom placed in a static 
electric field and illuminated by microwave fields, for the case 
with [panels (a)-(b)] and without [panels (c)-(d)] the parity 
symmetry II z = exp(i7rS z )P z compared with the predictions 
of random matrices ensembles. In panels (a) and (c) solid lines 
indicate the best fitting Izrailev distributions, while dashed 
and dotted lines correspond to GUE and GSE distributions 
respectively (in panel (c) the dotted line is hardly visible be- 
hind the solid one). Panels (b) and (d): solid and dotted 
(hardly visible behind the solid lines) lines correspond to nu- 
merical data and their best fits, while dashed lines indicate 
GUE and GSE predictions as indicated in the figure. 

Let us now consider the influence of the spin-orbit 
interaction in quantum case. The fine structure split- 
ting scales as n^ 3 . Thus, when collecting eigenvalues for 
different no manifolds (to improve the statistics) exter- 
nal fields have to be appropriately rescaled for the data 
to correspond to the same interesting physical situation 
in which external perturbations are comparable to the 
spin-orbit interaction. In effect we have diagonalized the 
quantum effective Hamiltonian for different hydrogenic 
manifolds in the range no = 50 — 59 keeping constant 
no-Bo = 10~ 6 , Fq/Eq and Fq 2 /E what ensures that the 
energy shift, with respect to the —l/2n 2 value, multiplied 
by n^ is of the same magnitude for each no manifold. 
For the case when H is invariant with respect to the par- 
ity symmetry, IT 2 we have collected single members from 
each Kramers pair in the range of "n 3 , x energy shift 
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FIG. 3. The best fitting Izrailev distribution parameter, 
for hydrogen atom in the static electric field and driven by 
the two microwave fields, as a function of the parameters of 
the fields. Panel (a): gradual breaking of the parity sym- 
metry, II 2 with a change of the amplitude Fq and for fixed: 
F$/E = 2, n E = 1(T 6 , 6 = tt/2 and <p = tt/4. For Fq = 
the symmetry is preserved, when Fq increases the symmetry 
is gradually broken and the statistics changes from the GUE 
to GSE type. Panel (b): switching off the spin-orbit inter- 
action - keeping the configuration of the external fields fixed 
(i.e. Fi/Eo = 2, Ff/Eo = 5, = 0.31tt and <p = tt/4) but 
increasing the fields amplitudes the spin-orbit interaction be- 
comes relatively weaker. This results in a change of the level 
statistics from GSE to GOE type. To calculate each point 
in the panels the effective Hamiltonian has been diagonalized 
for different hydrogenic manifold in the range no = 55 — 59. 



When the Hamiltonian is invariant with respect to the 
Il z transformation the nearest-neighbor spacing (NNS) 
distribution, P(s), obtained is fitted quite satisfactorily 
with the distribution appropriate for GUE see Fig. || 
top row. The solid line is a best fitting Izrailev distri- 
bution |l5|]l6f l with the parameter (3 = 1.91 (/3 gives the 
level repulsion exponent, i.e., P(s) cx s 13 for small spac- 
ing s) and with /N = 0.4, i.e. chi-squared divided by 
the number of levels N (there are about 14,000 levels in 
the data set). The spectral rigidities, A3, also reveals 
the behavior very close to GUE - the best fitting A3 
statistics corresponding to an independent superposition 
of Poisson and GUE spectra || results in the parameter 
value (a relative measure of the chaotic part of the phase 
space) q = 0.99. 

Breaking of the parity symmetry has a dramatic ef- 
fect on the quasienergies statistics as depicted in Fig. ^ 
- bottom row. This time both NNS distribution and A3 
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statistics obtained show a very good agreement with the 
predictions corresponding to GSE. In particular one ob- 
tains a quartic level repulsion, P(s) oc s for small spac- 
ings. More precisely we have got the fitted Izrailev dis- 
tribution repulsion parameter (3 = 3.86 with x 2 /N = 0.6 
and q = l from the best fitting A 3 statistics (now it cor- 
responds to an independent superposition of Poisson and 
GSE spectra |§). This time there are about 20,000 levels 
in the data set. 

Finally we would like to investigate gradual breaking of 
the parity symmetry Tl z as well as switching off the spin- 
orbit interaction J2|. To this end we have diagonalized 
the effective Hamiltonian changing the parameters of the 
system. In Fig. ^a we show how the best fitting Izrailev 
distribution repulsion parameter, /3 changes when Fq in- 
creases. That corresponds to a gradual breaking of the 
parity invariance - the transition from GUE type statis- 
tics to GSE one is clearly visible in the figure. Switching 
off of the spin-orbit interaction has been realized by in- 
creasing the microwave amplitudes which implies that the 
influence of the external fields becomes stronger while the 
spin-orbit interaction becomes relatively weaker. The re- 
sults, as shown in Fig. [|b, indicate the transition from 
GSE, through intermediate, to GOE type statistics. For 
negligible spin-orbit interaction additional unitary sym- 
metries are restored (e.g. IT = exp (in Si) where i = x,y 
with the property TI? = —1) which allows splitting the 
Floquet Hamiltonian into two identical real matrices 
10,0 an d consequently the GOE statistics is expected. 
Actually, without the spin-orbit interaction, the spin de- 
grees of freedom have no effect on a dynamics and can be 
eliminated. Then it becomes immediately clear that, be- 
cause of the time-reversal invariance, the system should 
reveal the GOE type statistics. 

To summarize, we have given an example of a real- 
istic physical system which may yield a quartic level 
repulsion and, more generally, have statistical proper- 
ties close to the Gaussian Symplectic Ensemble of RMT. 
This is possible by taking into account the spin-orbit in- 
teraction while using external perturbations (microwave 
and static) for creating a chaotic secular motion of the 
electronic ellipse. For the effect to be visible the spin- 
orbit and external perturbations have to be of compa- 
rable importance. This requires very weak external per- 
turbation since the spin-orbit coupling decreases rapidly 
with riQ. The resulting mean Floquet levels spacing 
within the n$ manifold is very small - on the edge of 
present experimental possibilities for unambiguous spec- 



troscopy. To improve experimental conditions one may 
employ hydrogen-like ions since the spin-orbit interaction 
increases with the charge of the nucleus and consequently 
mean level spacing could be greater. Importantly, the 
Rydberg states in question have large spontaneous emis- 
sion lifetimes so, at least in principle, the individual lines 
should be observable (and quartic level repulsion assures 
low probability of small spacings). 
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